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Abstract 

The proton analyzing power in carbon has been measured for energies of 82 to 217 MeV 
and proton scattering angles of 5 to 41°. The measurements were carried out using polarized 
protons from the elastic scattering ^H(e, reaction and the Focal Plane Polarimeter (FPP) 
in Hall A of Jefferson Lab. A new parameterization of the FPP p-C analyzing power was fit 
to the data, which is in good agreement with previous parameterizations and provides an 
extension to lower energies and larger angles. The main conclusions are that all polarimeters 
to date give consistent measurements of the carbon analyzing power, independently of the 
details of their construction and that measuring on a larger angular range significantly 
improves the polarimeter figure of merit at low energies. 
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1. Introduction 



In order to extract polarizations from asymmetries in proton-carbon scattering, accurate 
values of the analyzing power, A^, are needed. The analyzing power depends on the strength 



of the spin-orbit coupling and is a function of the proton kinetic energy, Tp, and of the ^o^ 
scattering angle, 9fpp. We have used a functional form similar to those used previously 
to obtain a new parameterization in the energy range of 82 to 217 MeV and proton scattering 



G3iar 



angles in carbon between 5 and 4 



previous parameterizations 



1°. The new parameterization is in good agreement with 



0, y, y, E 



5[ and provides an extension to lower energies and 



larger scattering angles. To date, the new parameterization has been used in three low 



energy polarization experiments 
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2. Method 

The experiment was carried out in Hall A of Jefferson Lab Q]. The primary elastic 
scattering ^B.{e,p) reaction provided the source of polarized protons for the calibration. 
The scattered proton polarization was determined independently of our knowledge of Ac. 

A continuous, polarized electron beam with polarization ranging from 80 to 85% was 
produced using a strained gallium-arsenide (GaAs) source 10|, The polarization in 



Hall A was limited to h = 38-41% due to multi-hall running. The beam helicity was flipped 
pseudo-randomly at 30 Hz, with negligible beam charge asymmetry between the two helicity 
states. The electron beam was accelerated to either 362 or 687 MeV and then scattered from 
a 15 cm long liquid hydrogen target. 

The scattered protons were detected in the left High Resolution Spectrometer (HRS) []J 
made up of one dipole and three quadrupole magnets. Vertical drift chambers were used to 
track the protons after deviation by the magnetic field of the dipole, allowing a determination 
of the momentum with a resolution of 2 x 10~^ [3]. The proton target quantities ytg 
(transverse position), dtg {dx/dz, where x is the displacement in dispersive plane and z is 
along the spectrometer axis), (ptg (dy/dz) and 6tg (fractional deviation of momentum from 
central trajectory) were reconstructed from the proton trajectory using the HRS optics 
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matrix. Triggering was performed using a coincidence between two scintillator planes, SI 
and S2, each made up of six panels. The scintillators SI and S2 are both located before the 
Focal Plane Polarimeter (FPP) carbon analyzer. 

The FPP, downstream of the VDCs and trigger panels, measured the recoil polarization of 
the protons through a secondary scattering of the protons from a carbon block "analyzer" . 
Spin-orbit coupling between the proton spin and angular momentum about the analyzer 
carbon nucleus leads to an asymmetry in the azimuthal scattering angle, (pfpp, reconstructed 
from the front and rear proton tracks. A detailed description of the polarimeter can be found 
in [l5|. The FPP measures the inclusive G{p,p') reaction at the low energies reported 
here, but also the inclusive C(p, X^*"), where X is a charged particle, for energies above pion 
production threshold. 

Cuts were made on the reconstructed proton target quantities of ytg < 65 mm, 6tg < 65 
mrad, (ptg < 38 mrad and 6tg < 0.045, in order to ensure that the protons originated in the 
target volume and traversed the HRS in a phase space region with high detection efficiency. 
To remove the majority of multiple scattering events due to Coulomb interactions between 
the proton and carbon nuclei, a cut on the FPP polar angle of 6fpp > 5° was made. A 
conetest cut was used to remove instrumental asymmetries 16|. In order to ensure that all 
FPP scattering events originated from within the carbon block, a cut was made around the 
location (along the spectrometer axis) of closest approach between the front and rear FPP 
tracks. A cut of 1 cm or less was also placed on the distance of closest approach between 
the front and rear tracks. 

Three carbon block thicknesses were used: 0.75, 2.25 and 3.75" . The carbon thickness was 
varied aeccdi,, .o GEANT fl Monte CaHo .tud.es cde. to ,na.i..e the «gu.e of .ent 
(FOM), which is a measure of how many events contribute to the scattering asymmetry 
and is given by the following integral over 6fpp or approximate summation over bins in 
Ofpp: 
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FOM = / e/pp(%pp)A,2(^/pp)rf%pp ^ Yl Ki^fpp^i (1) 

where Ac is the analyzing power and e/pp is the FPP efficiency, given by the fraction of events 
passing the target cuts and having tracks in the front FPP chambers which pass the FPP 
cuts. For the purposes of choosing carbon thicknesses, the Monte Carlo studies assumed Ac 



to follow the earlier parameterization of 



2j. A summary of the primary ^B.{e,p) reaction 



kinematics and FPP parameters can be found in Table [H 
3. Analysis 

The recoil transferred polarization was determined by a maximum likelihood method 
using the difference of the azimuthal distributions at the focal plane corresponding to the 
two beam helicity states. After spin transport through the spectrometer using COSY [isl . 
a differential algebra based code, the transferred polarization products at the target, hAcP'^ 
and hAcP'z-i were obtained. The ratio of these products was used to calculate the proton 
elastic form factor ratio, Ge^/Gmpi with the following equation jl^ : 



Gep ^ hAcP', Eq + Ec Oc 
Gm. hAcPl 2M, ''''2 



(2) 



where Eo and Eg are the energies of the incident and scattered electron, respectively, 6e is 
the electron scattering angle in the laboratory frame and Mp is the proton mass. Note that 
knowledge of the analyzing power and beam polarization are not required as they cancel. 



The form factor ratio calculated with Equation [2] was in turn used in t 
for the transferred polarization components at the target 19 
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le following derivations 
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22|: 



(3) 



, ^ E, + E, / . . tan^^ 

where r = Q'^/AM^, is the four-momentum transfer squared and e = [1 + 2(1 + 
r)tan^(Y)]~^ is the longitudinal polarization of the virtual photon. and P'^ calculated 
with equations [3] and HI respectively, were combined with the transferred polarization prod- 
ucts at the target obtained through azimuthal asymmetries, hAcP^ and hA^P'^^ to extract 



the analyzing power and the statistical uncertainty IJ], where h was taken to be the average 
beam helicity associated with the data. 

The extracted analyzing power data can be found in Tables [2] and [3] for incident electron 
beam energies of 362 and 687 MeV, respectively. The total systematic uncertainty due to 
uncertainties in beam polarization, spin transport, reconstruction of FPP scattering angles, 
proton momentum and beam energy is also shown. Table H] presents the weighted aver- 
age analyzing power, FPP efficiency and figure of merit, calculated using the summation 
approximation of equation [1] and the data in Tables [2] and [3l 

The functional form used to fit the extracted analyzing power is similar to the "low 
energy fit" used by McNaughton et al. E | . It is given by: 



where r = ppSm{dfpp) and pp is the proton momentum in GeV/c at the center of the carbon 
analyzer. Note that the functional form given satisfies the physical constraint that the 
analyzing power vanish at Ofpp = and 180°. The coefficients a, b, c and d are polynomials 
of the momentum. The last parameter, d, was added in order to improve the quality of fit. 
The coefficients are expanded as follows: 
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X = J2 MPp - Po)\ for X = a,b,c,d (6) 



i=0 



where Xi are the parameters of the fit. In order to obtain a numerically stable fit, the 
parameter po was set to the middle of the momentum range at the center of the carbon 
block, or 0.55 GeV/c. The quality of the obtained fits was largely insensitive to the choice 
of Pq, provided it was roughly in the middle of the momentum range being fit. The energy loss 
was approximated by dividing the carbon block into 1 cm layers and calculating the energy 



loss at each layer using a parameterization of the proton stopping power in carbon [23| . The 
procedure was repeated up to the center of the carbon block resulting in an average total 
energy loss. The entire data set found in Tables [2] and [3] was fit simultaneously and the 
parameters of this fit can be found in Table [5l 

To determine the statistical uncertainty of the fit the data points were shifted randomly 
within their statistical error bars and the resulting distribution was refit. The procedure 
was repeated 100 times and the maximum shift from central analyzing power value at each 
value of Tp and 6fpp was taken as the statistical uncertainty of the fit at that point. The 
systematic uncertainty of the fit was determined by fitting the distributions shifted both 
upward and downward by the systematic uncertainty of the data and the half width at 
each value of Tp and 6fpp was taken as the systematic uncertainty of the fit at that point. 
The statistical and systematic contributions to the total fit uncertainty at each point were 
then added in quadrature. The fit uncertainty was then parameterized as a function of Tp 
and 6fpp using the functional form of equation 5 and the parameters of the error fit can be 
found in Table [5l The uncertainty of the parameterized analyzing power can therefore be 
calculated for any kinetic energy at carbon center between 82 and 217 MeV and scattering 
angle between 5 and 41°. 
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4. Results 

Figures [U and [2] compare the new (solid), the "low energy" McNaughton (dashed), the 
"low energy" Aprile-Giboni [3| (dashed dotted), the Waters [4] (dotted) and the Pospischil [Sj 
(dashed double dotted) parameterizations to the new low energy data. Note that low energy 
data was also taken by M. leiri et al. 24| for kinetic energies at carbon center between 20 and 
84 MeV and scattering angles between 15 and 80° and by S.M. Bowyer et al. 25| for kinetic 
energies at carbon center between 120 and 200 MeV and scattering angles between 6 and 
23°. The fit uncertainty of the new parameterization, as described above, is denoted by the 
gray band. The "low energy" Aprile-Giboni parameterization was done for scattering angles 
of 5-20°, energies at carbon center of 90-386 MeV and a carbon block thickness of 3 cm. 
The Waters parameterization was done for scattering angles of 3.5-28°, energies at carbon 
center of 90-450 MeV and carbon block thicknesses of 3 and 6 cm. Note that Waters et al. 
quote a 10-20% uncertainty on their fit for energies between 90 and 100 MeV. The "low 
energy" McNaughton parameterization was fit to the world database at that time, including 
data from Aprile-Giboni et al. 3|] and Waters et al. ^, with carbon thicknesses ranging 
from 3 to 12.7 cm. It is valid for scattering angles between 5 and 20° and kinetic energies 
at the center of the carbon between 100 and 450 MeV. The parameterization of Pospischil 
et al., with a carbon block thickness of 7 cm, used a larger angular range of 15.7-43.1° but 
a narrow energy range of 160-200 MeV. It is shown only in the Tp = 171.1 MeV panel of 
Figure [21 

There is a good agreement between the new parameterization and the older ones in the 
energy/angle regimes for which they were intended, considering all fits were done for differ- 
ent polarimeters and for varying carbon block thicknesses. Consistency between different 
polarimeters has been observed before and is generally expected. This is in part because the 
energy range covered by most polarimeters in a single measurement is sufficiently small that 
the variation of analyzing power with energy is essentially linear. The new parameterization 
provides an extension down to an energy of 82 MeV and up to a scattering angle of 41°. 

The energy dependence of the weighted average analyzing power is shown for the angular 
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range of 5° < Ofpp < 20° in Figure [HI which again shows a very good agreement between 
the new data and the parameterizations in the energy/angle regimes for which the previous 
parameterizations were intended. 

A higher figure of merit results in more useful events and a lower statistical uncertainty on 
extracted polarization observables. The figure of merit as a function of maximum scattering 
angle is shown in Figures H] and [5] for beam energies 362 and 687 MeV, respectively. For 
the highest energies, the dependence of the figure of merit on maximum scattering angle is 
relatively flat beyond 25°. For the lowest energies, the figure of merit increases steadily as 
the maximum scattering angle is increased to 41°. By extending the maximum scattering 
angle of the parameterization, one obtains an improvement of up to 1.6 times in the figure 
of merit at low energies, as can be seen in Figure O 

5. Conclusion 

A new parameterization, similar to that of j^, was determined for p-C analyzing power 
in the 82 to 217 MeV energy region, and found to be in good agreement with the previous 
parameterizations in the energy and angular regimes for which they were intended. This 
corroborates that the various polarimeters built to date are measuring a property of inclusive 
proton-carbon scattering insensitive to details of how the polarimeters are constructed. The 
extension made for energies between 82 and 90 MeV and for scattering angles between 20 and 
41° increased the figure of merit by a factor of 1.6 at low energies. Any new polarimeters for 
low energy protons should take advantage of the increase in figure of merit from large-angle 
scattering in the polarimeter. 
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Ee 

(MeV) 


nlao 

(deg) 


Pp 

(MeV/c) 


T 
(MeV) 


Analyzer 
Thickness 
(inches) 


362 


28.3 


476.1 


82.2 


0.75 


362 


23.9 


503.2 


97.2 


0.75 


362 


18.8 


526.3 


105.3 


0.75 


362 


18.8 


526.3 


105.3 


2.25 


362 


14.0 


541.4 


118.7 


0.75 


687 


47.0 


583.9 


126.9 


2.25 


687 


44.2 


621.1 


134.9 


3.75 


687 


40.0 


674.0 


171.1 


3.75 


687 


34.4 


740.3 


216.6 


3.75 



Tabic 1: Kinematics of the primary ^H(e, p) reaction and FPP parameters. E^, Oiab and Pp are the beam 
energy, proton lab angle and proton central spectrometer momentum, respectively. Tp is the proton kinetic 
energy at the center of the carbon analyzer. Note that two carbon thicknesses (0.75 and 2.25") were used 
for the 6*'°^ = 18.8° kinematic setting. 
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Figure 1: Analyzing power data and parameterization for beam energy E^, = 362 MeV. Error bars shown 
are statistical and systematic uncertainties added in quadrature. The dashed line denotes the "low energy" 
McNaughton parameterization 0, the dashed dotted the "low energy" Aprile-Giboni parameterization Q , 
the dotted line the Waters parameterization Q , the solid line the new parameterization from this work and 
the gray area the error band. 
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Tp 6 fpp y4c,i± Stat ± syst efpp^i 

(MeV) (deg) (%)_ 

82.2 



97.2 



105.3 



118.7 



7 


0.143 


± 


0, 


.010 


± 


0, 


.007 


0.33 


11 


0.177 


± 


0, 


.008 


± 


0, 


.008 


0.42 


15 


0.180 


± 


0, 


.008 


± 


0, 


.008 


0.45 


19 


0.190 


± 


0, 


.009 


± 


0, 


.009 


0.35 


23 


0.186 


± 


0, 


.011 


± 


0, 


.008 


0.22 


27 


0.206 


± 


0, 


.015 


± 


0, 


.010 


0.12 


31 


0.261 


± 


0, 


.020 


± 


0, 


.013 


0.07 


35 


0.320 


± 


0, 


.026 


± 


0, 


.016 


0.04 


39 


0.371 


± 


0, 


.031 


± 


0, 


.017 


0.03 


7 


0.212 


± 


0, 


.010 


± 


0, 


.009 


0.28 


11 


0.234 


± 


0, 


.007 


± 


0, 


.010 


0.41 


15 


0.245 


± 


0, 


.007 


± 


0, 


.010 


0.41 


19 


0.268 


± 


0, 


.009 


± 


0, 


.012 


0.29 


23 


0.283 


± 


0, 


.011 


± 


0, 


.012 


0.17 


27 


0.252 


± 


0, 


.015 


± 


0, 


.011 


0.10 


31 


0.216 


± 


0, 


.019 


± 


0, 


.011 


0.06 


35 


0.208 


± 


0, 


.023 


± 


0, 


.009 


0.04 


39 


0.192 


± 


0, 


.033 


± 


0, 


.010 


0.03 


7 


0.124 


± 


0, 


.008 


± 


0, 


.009 


1.93 


11 


0.258 


± 


0, 


.009 


± 


0, 


.017 


0.68 


15 


0.290 


± 


0, 


.009 


± 


0, 


.018 


0.58 


19 


0.319 


± 


0, 


.011 


± 


0, 


.020 


0.37 


23 


0.328 


± 


0, 


.014 


± 


0, 


.020 


0.20 


27 


0.266 


± 


0, 


.017 


± 


0, 


.017 


0.11 


31 


0.228 


± 


0, 


.021 


± 


0, 


.014 


0.06 


35 


0.179 


± 


0, 


.024 


± 


0, 


.010 


0.04 


39 


0.080 


± 


0, 


.032 


± 


0, 


.007 


0.03 


7 


0.249 


± 


0, 


.008 


± 


0, 


.012 


0.26 


11 


0.318 


± 


0, 


.007 


± 


0, 


.014 


0.39 


15 


0.361 


± 


0, 


.007 


± 


0, 


.016 


0.36 


19 


0.379 


± 


0, 


.009 


± 


0, 


.017 


0.25 


23 


0.360 


± 


0, 


.011 


± 


0, 


.016 


0.16 


27 


0.270 


± 


0, 


.013 


± 


0, 


.014 


0.11 


31 


0.184 


± 


0, 


.015 


± 


0, 


.008 


0.09 


35 


0.140 


± 


0, 


.017 


± 


0, 


.009 


0.06 


39 


0.094 


± 


0, 


.023 


± 


0, 


.008 


0.03 



Table 2: Extracted analyzing power for electron beam energy ~ 362 MeV. 
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Tp 6 Ac^,± Stat ± syst ejpp^i 

(MeV) (deg) (%^ 

126.9 



134.9 



171.1 



216.6 



7 


0.244 


± 


0, 


.006 


± 


0.016 


1.34 


11 


0.379 


± 


0, 


.007 


± 


0.024 


1.24 


15 


0.433 


± 


0, 


.008 


± 


0.027 


0.97 


19 


0.453 


± 


0, 


.010 


± 


0.029 


0.59 


23 


0.415 


± 


0, 


.012 


± 


0.026 


0.34 


27 


0.339 


± 


0, 


.016 


± 


0.022 


0.21 


31 


0.270 


± 


0, 


.018 


± 


0.017 


0.15 


35 


0.198 


± 


0, 


.023 


± 


0.013 


0.10 


39 


0.120 


± 


0, 


.029 


± 


0.008 


0.07 


7 


0.179 


± 


0, 


.004 


± 


0.008 


2.35 


11 


0.403 


± 


0, 


.005 


± 


0.016 


1.47 


15 


0.443 


± 


0, 


.005 


± 


0.017 


1.14 


19 


0.480 


± 


0, 


.007 


± 


0.019 


0.68 


23 


0.436 


± 


0, 


.009 


± 


0.018 


0.38 


27 


0.353 


± 


0, 


.011 


± 


0.014 


0.21 


31 


0.295 


± 


0, 


.016 


± 


0.011 


0.13 


35 


0.236 


± 


0, 


.024 


± 


0.009 


0.08 


39 


0.181 


± 


0, 


.022 


± 


0.007 


0.05 


7 


0.394 


± 


0, 


.006 


± 


0.013 


1.38 


11 


0.559 


± 


0, 


.006 


± 


0.019 


1.28 


15 


0.600 


± 


0, 


.007 


± 


0.020 


0.93 


19 


0.493 


± 


0, 


.009 


± 


0.018 


0.60 


23 


0.334 


± 


0, 


.010 


± 


0.013 


0.44 


27 


0.203 


± 


0, 


.011 


± 


0.007 


0.33 


31 


0.109 


± 


0, 


.012 


± 


0.004 


0.25 


35 


0.048 


± 


0, 


.014 


± 


0.002 


0.18 


39 


-0.043 ± 0.036 ± 0.004 


0.11 


7 


0.537 ± 0.006 ± 0.030 


1.23 


11 


0.677 ± 0.006 ± 0.037 


1.15 


15 


0.584 ± 0.007 ± 0.033 


0.85 


19 


0.380 ± 0.008 


± 0.022 


0.67 


23 


0.219 ± 0.008 


± 0.012 


0.59 


27 


0.134 ± 0.008 


± 0.008 


0.53 


31 


0.047 ± 0.009 ± 0.003 


0.41 


35 


0.045 ± 0.048 


± 0.017 


0.30 


39 


-0.060 ± 0.021 ± 0.003 


0.18 



Table 3: Extracted analyzing power for electron beam energy E,, ~ 687 MeV. 
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T 


^fpp 


< Ac > ±stat ± syst 


FOM 


(MeV) 


(%) 




(%) 


82.2 


2.03 


0.186 ± 0.010 ± 0.009 


0.07 


97.2 


1.80 


0.243 ± 0.010 ± 0.010 


0.11 


105.3 


4.01 


0.239 ± 0.011 ± 0.015 


0.20 


118.7 


1.72 


0.307 ± 0.009 ± 0.014 


0.17 


126.9 


5.02 


0.354 ± 0.009 ± 0.023 


0.66 


134.9 


6.50 


0.336 ± 0.006 ± 0.014 


0.81 


171.1 


5.49 


0.423 ± 0.009 ± 0.015 


1.16 


216.6 


5.91 


0.391 ± 0.010 ± 0.023 


1.31 



Table 4: Weighted average analyzing power (< Ac >), FPP efficiency (e/pp) and figure of merit (FOM). 
The efpp values are sums of e^^^ over all bins in 9 fpp in Tables [5] and [31 The FOM values were calculated 
using the summation approximation of equation [T] and, along with < Ac >, the entire angular range of 
5° < 0fpp < 41°. 



Coefficient 


Value 


Unit 




Po 


0.55 


(GeV/c) 




ao 


4.0441 


(GeV/c) 


-1 


ai 


19.313 


(GeV/c) 


-2 


a2 


119.27 


(GeV/c) 


-3 


as 


439.75 


(GeV/c) 


-4 


a.4 


9644.7 


(GeV/c) 


-5 


bo 


6.4212 


(GeV/c) 


-2 


bi 


111.99 


(GeV/c) 


-3 


b2 


-5847.9 


(GeV/c) 


-4 


bs 


-21750 


(GeV/c) 


-5 


b4 


973130 


(GeV/c) 


-6 


Co 


42.741 


(GeV/c) 


-4 


Cl 


-8639.4 


(GeV/c) 


-5 


C2 


87129 


(GeV/c) 


-6 


cs 


8.1359 xlO^ 


(GeV/c) 


-7 


C4 


-2.1720 xlO^ 


(GeV/c) 


-8 


do 


5826.0 


(GeV/c) 


-6 


di 


2.4701x10^ 


(GeV/c) 


-7 


d2 


3.3768 xlO^ 


(GeV/c) 


-8 


ds 


-1.1201x10^ 


(GeV/c) 


-9 


d4 


-1.9356x10^ 


(GeV/c) 


-10 



Table 5: Parameters of analyzing power polynomial fit; reduced of the fit is 1.39 with a of 71.0 and 
51 degrees of freedom. 
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Error 
Coefficient 


Value 


Unit 




Po 


0.55 


(GeV/c) 






0.34687 


(GeV/c) 


—1 


1 


3.8266 


(GeV/c) 


—2 




11.635 


(GeV/c) 


—3 


ag 


-314.18 


(GeV/c) 


—4 


3-4 


1015.1 


(GeV/c) 


—5 




0.21579 


(GeV/c) 


—2 




288.29 


(GeV/c) 


—3 


bi 


-13736. 


(GeV/c) 


—4 


b^ 


-1.5821 xlO^ 


(GeV/c) 


— 


b2 


2.2523 xlO'^ 


(GeV/c) 


-6 


ce 


938.40 


(GeV/c) 


-4 




-4961.8 


(GeV/c) 


-5 




4.0190 xlO^ 


(GeV/c) 


-6 




3.3034 xlO^ 


(GeV/c) 


-7 


cl 


-5.8401x10^ 


(GeV/c) 


-8 




-4910.9 


(GeV/c) 


-6 




1.5133 xlO^ 


(GeV/c) 


-7 




-1.3127 xlO^ 


(GeV/c) 


-8 




-2.7446x10^ 


(GeV/c) 


-9 


U4 


2.8556x10® 


(GeV/c) 


-10 



Table 6: Parameters of fit to the uncertainty in analyzing power parameterization. 
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Figure 2: Analyzing power for beam energy = 687 MeV. Error bars shown are statistical and systematic 
uncertainties added in quadrature. The dashed line denotes the "low energy" McNaughton parameteriza- 
tion 0], the dashed dotted the "low energy" Aprile-Giboni parameterization 0], the dotted hne the Waters 
parameterization the dashed double dotted Hne the Pospischil parameterization the sohd hne the 
new parameterization from this work and the gray area the error band. 



o 




100 150 200 

Tp (MeV) 

Figure 3: Weighted average analyzing power as a function of proton energy at the center of the carbon (Tp) 
for 5° < dfpp < 20°. Error bars shown are statistical and systematic uncertainties added in quadrature. 
The dashed line denotes the "low energy" McNaughton parameterization the dashed dotted the "low 
energy" Aprile-Giboni parameterization 0], the dotted line the Waters parameterization 01, the solid line 
the new parameterization from this work and the gray area the error band. Note that oscillations in the 
data near Tp — 130 MeV arc due to changes in the analyzer thickness. 
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Figure 4: Figure of merit as a function of maximum scattering angle for beam energy Eo = 362 MeV. Solid 
line denotes proton kinetic energy at center of carbon analyzer Tp = 82.2 MeV, dashed line Tp = 97.2 MeV, 
dotted line Tp = 105.3 MeV and dashed dotted hue Tp = 118.7 MeV. Note that the figure of merit for Tp = 
105.3 MeV (0.75 and 2.25" carbon blocks) is larger than for Tp = 118.7 MeV (0.75" carbon block) because 
at these energies the figure of merit is greater for a 2.25" analyzer thickness than a 0.75" one. 
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Figure 5: Figure of merit as a function of maximum scattering angle for beam energy Eq = 687 MeV. Solid 
line denotes proton kinetic energy at center of carbon analyzer Tp = 126.9 MeV, dashed line Tp = 134.9 
MeV, dotted hue Tp = 171.1 McV and dashed dotted hue Tp = 216.6 MeV. 
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Figure 6: Figure of merit ratio for two 9fpp ranges {FOM{h — 41°) / FOM{b — 20°)) as a function of proton 
energy at carbon center (Tp). Uncertainties were calculated using the analyzing power error parameterization 
and the coefficients in Table [S] 
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